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We obtain through a Matrix Product Ansatz (MPA) the exact solution of the most general N- 
state spin chain with U(1) N symmetry and nearest neighbour interaction. In the case N = 6 this 
model contain as a special case the integrable SO(6) spin chain related to the one loop mixing 
matrix for anomalous dimensions in N — 4 SYM, dual to type II B string theory in the generalised 
Lunin-Maldacena backgrounds. This MPA is construct by a map between scalar fields and abstract 
operators that satisfy an appropriate associative algebra. We analyses the Yang-Baxter equation 
in the N = 3 sector and the consistence of the algebraic relations among the matrices defining the 
MPA and find a new class of exactly integrable model unknown up to now. 
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INTRODUCTION 

Since the discovery of the relation between the planar dilatation operator of the J\f = 4 Super Yang-Mills with 
an integrable so(6) quantum spin chains [I], integrability has played a prominent role in the exploration of the 
Maldacena's correspondence between the JV = 4 Super Yang-Mills with IIB string theory in AdS^ x S 5 spaces 
[21 El H]. The study of the planar dilatation operator's integrability is very important because it not only enable us 
to test the Maldacena's correspondence as it is an generator of nontrivial integrable models. Exactly solvable models 
are of interest in both physics and mathematics since the pioneering work of Hans Bethe [5]. The Bethe ansatz 
and its generalisations emerged over the years as a quite efficient and powerful tool for the exact solution of the 
eigcnspectrum of a great variety of one dimensional quantum chains and two-dimensional transfer matrices (see, e.g., 
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[El HI IS] for reviews) . According to this ansatz the amplitudes of the eigenfunction are expressed by a nonlinear 
combination of properly defined plane waves. On the other hand, in the last two decades several different ansatz 
were introduced in the literature under the general name of matrix product ansatz (MPA). The first formulation 
was done for the description of the ground-state eigenfunction of some special non-integrable quantum chains, the so 
called valence-bond solid models [TUl EJ EH H3] ■ The MPA becomes also a successful tool for the exact calculation 
of the stationary probability distribution of some stochastic one dimensional systems |14l 1151 116] . An extension of 
this last MPA, called dynamical MPA was introduced in [17l [18] and extended in [19 . This last ansatz gives the 
time-dependent probability distribution for some exact integrable systems. The MPA we are going to use in this 
paper, in order to solve the U(1) N quantum spin chain, was introduced in [20] HB 1221 123] ■ This ansatz was applied 
with success in the evaluation of the spectra of several integrable quantum Hamiltonians [20] [21] [22] , transfer matrices 
[23] [551125] and the time-evolution operator of stochastic systems [23] . According to this ansatz, the amplitudes of the 
eigenfunctions are given in terms of a product of matrices where the matrices obey appropriated algebraic relations. 
In the case of the Bethe ansatz the spectral parameters and the amplitudes of the plane waves are fixed, apart from a 
normalisation constant, by the eigenvalue equation of the Hamiltonian or transfer matrix. On the other hand, in the 
MPA the eigenvalue equation fixes the commutation relations of the matrices defining the ansatz. In such case the 
spectrum of the Hamiltonian or transfer matrix, and the corresponding eigenfunctions, can be computed in a purely 
algebraic way. 

In the present paper we study the most general quantum spin chain with U(1) N symmetry with nearest neighbour 
interaction thought a MPA in the generalised Lunin-Maldacena backgrounds (in this case we are only interested in 
zero momentum eigenstates, the more general case, including all momentum states, will be presented in |27j). This 
model contains as a special case the integrable so(6) spin chain related to the one loop planar dilatation operator of 
the Af = 4 Super Yang-Mills [1J. We analyse the Yang-Baxter equation in the N = 3 sector and the consistence of 
the algebraic relations among the matrices defining the ansatz and find that the solutions are separated in two class. 
In the first (class A) we obtain the models presented in [2E1 1201 1301 IS] , m the context of one loop dilatation operator, 
as well as in condensed matter physics and stochastic models [111 HU H31 1321 1331 131] • In the second class (class B) we 
obtain the model presented in |19j for the stochastic problem of fully asymmetric diffusion of two kinds of particles. 
In this last sector we also find a new type of integrable model unknown up to now. Our solution generalises the 



previous results obtained in (14 1 through the coordinate Bethe ansatz for U(l) quantum spin chains. The analyses 
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of the model for the full sector with U(1) N symmetry will be presented elsewhere [27] . 

THE U(1) N SPIN CHAIN IN THE GENERALISED LUNIN-MALDACENA BACKGROUNDS 

The AdS/CFT conjecture relates operators, states, correlation functions and dynamics between the J\f = 4 Super 
Yang- Mills with IIB string theory in AdS^ x S 5 spaces. One of the most important results of this conjecture predicts 
that the spectrum of scaling dimension D of gauge invariant operators, in the conformal field theory, should coincide 
with the spectrum of energies E of string states. Furthermore, this correspondence relates the weak coupling constant 
regime, in the gauge theory, with the strong coupling constant ones, in the string theory. Currently it is only possible 
to test the Maldacena's conjecture for a limited class of operators, and only in the case of large A values of the 't 
Hooft coupling constant. Operators where this conjecture can be tested are called BMN (Berenstein, Maldacena e 
Nastase) operators 35J. These operators are local primary operators in S 5 with large J charge, in the planar limit 
(large A). These operator are given by the trace of a product of a large number the six scalar fields of the theory. The 
most general operator O in this class is given by 



where oij — 1, N, with N — 6 and j = 1, ...L, and L >> 1 is the number of fields. Instead of restrict our study 
to the case N — 6, we will consider in this work N general. The planar dilatation operator D, that is one of the 
generators of the conformal algebra and it give us the scaling dimensions of local operators, and consequently give us 
the energies of the string's states, is given by: 



where g = s ™ 2 , and the operator L counts the number of scalar fields on the operators O |l| . The operator H is 
the one loop planar dilatation operator. This operator H was identified [T| with an exactly integrable uni dimensional 
SO(6) quantum spin chain Hamiltonian. By this expansion in loops ^ it can perturbativelly compare the string 
theory in the BMN limit to the field theory. On the other hand, this equivalence between the dilatation operator with 
integrable spin chain makes possible to use the powerful Bcthc ansatz to diagonalise the one- loop dilatation operator 
[I]. By diagonalising this operator through the Bethe ansatz, it is obtained the scaling dimensions of any state of 
local operators in the conformal field theory in a purely algebraic way. The integrability also guarantees that these 
scaling dimensions can be exactly obtained for all gauge group, and in particular, for the U(N) group [36j . On the 
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side of string theory, the integrable structures was discovered by the observation that Green Schwartz superstrings in 
AdS§ x S 5 has an infinite set of conserved non-local charges [37] . These results have made possible the comparison 
between gauge theories and string theories in the plane wave limit |38j . furthermore, the integrability enable us to 
make very precise tests for the correspondence AdS/CFT. 

The most general one-loop dilatation operator is related with an integrable spin model with long-range couplings 
[31?] . In the particular case of only two types of scalar fields on the BMN operators ^ (case N = 2), this integrable 
Hamiltonian reduces to the famous XXX model pQ. On the other hand, the case of N — 3 was studied by Frolov e 
Tseytlin [40] and it corresponds to a spin-1 spin chain. The most general spin-1 model in condensed matter physics 
[2"T] was obtained by us in a recent work, but it is not studied with the objective of classifying all possible solutions. 
In the context of the one-loop dilatation operator the solution, on class A presented here, was obtained in [28 . 

We consider here the most general iV-state spin chain with nearest neighbour interaction, periodic boundary con- 
dition, zero momentum eigenstates [46] , and U(1) N symmetry. The N possible states configurations of a given site is 
related to the N different types of scalar fields in |lj. The U(1) N symmetry imply that the Hamiltonian describing 
the time evolution of this spin chain conserves the number of states of each type. By denoting the basis of states at 

a given site as |a) (a = 1, ...,N), the Hamiltonian in a periodic lattice with L sites takes the form 

l I n n \ 

^ = E E E KiEt-tfA ), (3) 

3=1 \a^/3=l a, 0=1 j 

where E a @ are N x N Wyel matrices with elements (i? Q P) . = 5 a ,*^8j (a, (3=1, N). While the first term in the 
right hand side of (jij) acts over neighbour sites exchanging its configuration |a) ® \(3) — ► \j3) ® \a) with rate ^, the 
second one is a diagonal operator with weight ^. The eigenfunction for (jij) can be construct as 

\0(*,N))= Yl *%- c "'\ai,:.,aL) (ai = l,...,Aa (4) 

ai, ...,(*£ 

where the amplitudes \fy<^ - aL are related to the amplitudes in (JT|), for the operator O, and the symbol (*) in the sum 
denotes the restriction to the sets {a\, a^} with the same number n a of spins in configuration a. 



THE MPA 

In order to formulate a MPA for the Hamiltonian ^ , we make a one-to-one correspondence between configurations 
of spins (or scalar fields in the SYM) and products of abstract matrices. This matrix product is construct by 
making a correspondence between sites with spin configuration a = 1,..,,N and a matrix A^ a \ Our MPA asserts 



that the components of the amplitude of the eigenfunction \jj c ^ - aL m Q are obtained by associating them to a 
products of these matrices A^ a \ Actually (a — 1, N) are abstract operators with an associative product. A 
well defined eigenfunction is obtained, apart from a normalisation factor, if all the amplitudes arc related uniquely. 
Equivalcntly, in the subset of words (products of matrices) in the algebra containing n a matrices (a = l,...,N, 
ni + • • ■ rijv = L) there exists only a single independent word. The relation between any two words gives the ratio 
between the corresponding amplitudes of the components of the eigenfunction |0(^, N)). To formulated the ansatz 
we can choose any uniform operation on the matrix products that gives a non-zero scalar. For a quantum chains with 
periodic boundary conditions the trace operation is a convenient chose to produce this scalar |47j . The amplitudes 
$^- Qi in (|4) takes the form: 

$J'"' 01 =T r {A^A^ ■■■A^) ( aj = l,...,N). (5) 

It is obvious that the N states \a...a) (a — 1, ...,N) are all eigenstates of the Hamiltonian (|3j. In the following we 
shall choose |l...l) as our reference state. The Hamiltonian (J3J) when applied to the components of the eigenfunction 
Q where we do not have spins configurations \a) (a — 2, ...,N) at nearest neighbour sites give us the constraints, for 
the amplitudes (jsj, 

e n Tr(A Xl ~ 1 A^ ai ^ A X2 ~ Xl ~ 1 ■ ■ ■ A^ aj -^ A X3 ' X3 - 1 ~ 1 A^ a ^ A X3+1 ~ Xj ~ 1 A^ aj+l) ■ ■ ■ A^^ A L ~ Xn ) 

n 
3=1 

QLj 1 V / 

n 

+[(L - 2n)T\ \ + (Kl I + T l 2)} r Tr(A Xl - 1 A {ai) A X2 - x ^ 1 ■ ■ ■ A {a " } A L ~ x ") (aj = 2, N), (6) 



1=1 



where e n is the energy of the eigenfunction A = A^\ n = ri2 + "3 + ■ ■ ■ + njsr, and the position in 

the spin chain where we have a state configuration \a =/= 1). A convenient solution of this last equation is obtained by 
identifying the matrices A^ (a — 2, ...,N) as composed by spectral-parameter-dependent matrices [201 [2TJ [22] . The 
distinguibility of states configurations allows two types of solutions. The standard solution (class A) is obtained if 
each of the matrices A^ (a = 2, TV) is composed of n = na +113 + • • • + ft.iv spectral parameter dependent matrices 
[2T1 125] . The second class of solutions (class B) is obtained if matrices A^ (a = 2, ...,N) with different a value are 
composed of by distinct sets of spectral parameters matrices [2T] , 
Solutions of class A 
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In this case, the matrices (a — 2, N) can be written in terms of the matrix A and n = 712 + TI3 + • • • + tin 



spectral parameter dependent matrices Ai. |48j : 



^ (Q) =E4: } A (a = 2,...,N), 



(7) 



where the matrices A.jj.. (a — 2, ...,N) satisfy the following commutation relations with the matrix A: 



A ( g ) A = g(a)e^AA%>, (a=2,...,N), (j=l,...,n), 



(8) 



the parameters fcj (j = 1, n) are in general complex numbers unknown a priori, and <?(a) is a constant. The energy 
e„ is obtained by inserting ([7]) in by using ([8| and imposing that s n is a symmetric function on the spectral 
parameters: 

n N 



2^ ikj 



E n « (ri2 + r^) + (^- 2 «) r i^ 



(9) 



Q = 2 



where we need to impose 



m = S| (3(2) = 1, a = 3, ...,7V). 



(10) 



pa 1 pi z 
1 1 a 1 2 1 

The relations coming from the eigenvalue equation for configurations where we have two spins configurations \a 
(a = 2, ...,N) at nearest neighbour sites are given by 



n 

E [ r 2i + r? 2 < 



ma I 1 pi a pi I pa a^ ifc. 



3,/=2 



4 a) 4?^o, 



(ii) 



E [rj? + r? 2 ^ fe > + ^ + CTSi + rJg-rii-rs J)e<*] = E 7^ r " (« * 0, 



where we have used ( 5 ) , The relations Jill fix the algebraic relations among the matrices A^ (a — 2, N) 

N 



47< ] = E S^MA^Af^ A^A^ = (l^j=l,..,n), 

a',/3'=2 

where the nonzero structural constants S% a ,(kj,k[) are given by: 



(12) 



Sa a (%i 

a(%' fa) 



pi 2 , p2 l„i(kj+ki) I /™ 1 i fl a _ fH _ pa a\ ifc; 
1 21+ 1 12 e t l 1 ttl tl la Ml 1 aa^ 

pi 2 I p2 lpi(kj+ki) i (rn 1 j. fl a _ rl 1 _ r Q »1A ' 



( aj L(3 = 2, ...,N) 



(13) 



c ai/9 (* is fez) - ri 2 + r 2 Je'^+fcO + (r« } + r* g - r\ \ - r a a g)e 4 *' . 



Relations (|8j) and (12 1 define completely the algebra whose structural constants are the S-matrix given by (13). It 
is important to mention that in the sector N — 2 (only two types of scalar fields (f>\ and (f>2, for example) the 
Hamiltonian (|3| reduces to the well known asymmetric XXZ model |41j . In this particular case the 5-matrix is a 



diagonal matrix were Sjj(kj,ki) in ( 13 1 is the only nonzero clement. As a consequence, the algebra defined by ( 12 1 is 



associative for TV = 2. For N general the S-matrix is not diagonal and the algebra ( 12 1 is not associative for arbitrary 



values of ^ and ^. For arbitrary amplitudes we have in our matrix product ansatz (|5|) a product of n matrices 



A^\ Our ansatz will be valid only if the relations (|12| provide a unique relation among these products, otherwise 



the eigenfunction (|4| is not properly defined. This means, for example, that the products • • ■ A^A^jA^ ■ ■ ■ and 
• • • A^A^A^ ■ ■ ■ should be uniquely related. Since we can relate then either by performing the commutations in 
the order a/37 ~~ > P a l ~ > Pl a ~ * 7 P a or a Pl ~ > a lft ~ > "fctfi — > 7/3a, the structure constants 5"y of the algebraic 



relations (12 1 should satisfy 



N N 



j2 s^:(k 1 ,k 2 )sj^(k l> h)sj:}:,(k 2 ,h)= ^ s^^ik^s^'^Msy^M, (wj 

7,7 / ,7 // — 2 7, 7', 7"— 2 



for a, a', a", /3, /3', (3 — 2, 6. This last constraint is just the Yang-Baxter relation [6, 42J of the 5-matrix defined in 
( j 1 3 j ) . Actually the condition |l4| is enough to ensure that any matrix product of spectral matrices {A^"'} is uniquely 
related and it implies the associativity of the algebra of the matrices {^t^} ■ In this case the relations coming from the 
eigenvalue equation for configurations where we have more than two spins configurations \a) (a = 2, N) at nearest 



neighbour sites are automatically satisfied by the algebraic relations (jSl) and (12 1 among the matrices defining the 



ansatz, no new constraints appear from these relations. On the other hand, the Yang-Baxter relation (14 1 produces 



strong constraints in the allowed couplings ^ and ^ of the Hamiltonian (j3J). In the sector N = 3 (tree types of 
scalar fields tfi±, 4>2 and (f>3, for example) the Hamiltonian ^ reduces to an asymmetric spin-1 quantum Hamiltonian 



with U(l) symmetry. The spin-1 model with g(2) = g(3) = 1 in (10 1 and U(l) symmetry was obtained by us in a 



recent work [5T] in the context of condensed matter physics, but it is not studied with the objective of classifying all 



possible solutions of the Yang-Baxter equation (14). In the generalised Lunin-Maldacena backgrounds, the solutions 



of ( 14 1 was classified in |28j for hermitian Hamiltonians (|3| . In the present paper we obtain the most general solution 



of (14). These solutions generalises the previous results obtained in [51] [5S] and can be obtained by a systematic 



investigation with the aid of Maple by purely analytical means. We find the following types of solutions: 



A.i) r„"r"^ — r^rf^ — r^r^ f — t a pitp a x, t aal — tppi — o, 

A ry\ pl otiq 1 pl 0-pP 1 -p/3 a-pa , , , , p. 

Sl-^) 1 Ql 1 la~ i ^l i l/3 — 1 a J3 1 (J a ~ <<a/31 t/3al , ^al — ^aa/3 — U, 

A.3) r„"r" * = r^rf^ = r^r^ f = t a pitp a x, hp a = taip = tp a i, 

aa) r^"r"^ = rl^rf 1 r^o = r2f = o, t a 0i = t 0a i = t i3i = t aa i, 



as) r^ = s rfi, 
A6) r^ = s rf£, 
a?) ri? = fl r^, 



pi a pi /3 r> pQ /3 . . 

1 a 1 — /31 — ' a — ' L a01 — <<aal) 



pl a _ pi _ n pQ ri j. _ j 

1 a 1 — /31 — ' /3a — ' t '/3al — ^aal) 



(15) 



1 1 ■ 



t-ii a -ril [3 T-ra i -rip i ri a p r\ , 



^0 1 
- 1 P 



^P 1 



o, r~ f — o, t 0( gi — f^i, 



a /3 



4 n\ pl a pa 1 pl P pP 

-^■ a J L al~ L la— L 01— L l 



^P 1 



0. 



where a, f3 = 2, 3 with a ^ /3 and t a , ( g, 7 = ^' + ^ - T%, } - (a', /?', 7 = 1,2,3). The first four models 

A.1-A.4 and A. 9 contain as a special case (for hermitian hamiltonians ([3])) the solutions obtained in [251 HHJ GUI EI] , m 
the context of one loop dilatation operator, as well as in condensed matter physics and stochastic models [THJ [211 1231 
32, 33, 34J. The models A.5-A.8 generalises the model presented in [19] for the stochastic problem of fully asymmetric 
diffusion of two kinds of particles. 

In order to complete our solutions through the MPA ^ we should fix the spectral parameters, or momenta, 
k\, . . . , k n . Using the algebraic relations (7l and (8 1 an arbitrary amplitude is proportional to Tr A^ 1 ' • • • A^"^ A L 



W 7 _ p7 0' 



The cyclic property of the trace and the commutation relations (|8j) and ( 12 1 give us 

N 



Tr 



.. . A^ n) A L 



ik-j L 



(ai,...,a n \T^ K,..., a'„)Tr 



■ ■ ■ n yu: ™> F 



(16) 



a , ... ,a' —2 



where we have used the identity (see (fl3"|) 



and 



E 

at'J.a" 



(kj, kj) 



(17) 



a",...,a^ 

where 4>(a'{) — g(a") L , is a (A — 1)" x (N — 1) "-dimensional transfer matrix of an inhomogeneous vertex model 



(inhomogeneities {fc;}) with Boltzmann weights given by ( 13 1 . The model is defined on a cylinder of perimeter n with 



a seam along its axis producing the twisted boundary condition 



S%& 1 (k n) k s ) = S%%(k n ,k i Wc®. 



(19) 



Finally relation (16) with (19 1 give us the constraints for the spectral parameters: 



-ikrj L 



= A^ (i = l,...,n), 



(20) 



where A^(kj,{h}) are the eigenvalues of the transfer matrix (1181. The condition (20 leads to the problem of 



evaluation the eigenvalues of the inhomogeneous transfer matrix (18). This can be done through the algebraic Bcthc 



ansatz [43 or the coordinate Bethe ansatz (see [44] and [45] for example). 
Solutions of class B 

In this case, the matrices (a = 2, ...,N) with different a value are composed of by distinct sets of spectral 
parameters matrices. For simplicity we are going to consider that there are two sets of spectral parameters. For a 



given (5 = 2. N, we have n a = n 2 + ■ ■ ■ + n$ parameters k^ (j = 1, n a ) and rib = n«5+i + ■ ■ 
^4 (j = 1, ni,). In this case the matrices can be written as: 

n a rib 

A(a) = Y. A % A ' («=2,~,o, ^ (Q) = E^S> A ' («=5+i,...,iv), 



njv parameters 



(21) 



3=1 



3 = 1 



where the matrices A^ 1 {a — 2, ...,N) satisfy the following commutation relations with the matrix A: 



A%A = g(a)e**rAA™, (a = 2, A), 



■ik 



(a) 



(22) 



with fcj a ' = k^ for a < 5 and k^ — k^ for a > 5. The energy e n is obtained by inserting (21 1 in (jGJ), by using (22 1 
and imposing that e n is a symmetric function on each sets of spectral parameters: 



e n = £ (r? \e*r + \e- ik r) + £ (rf£Je tt J iJ + rj*+i« 



JV 



X>«( r iS+rsl) + (£-2n)rii, (23) 



where we need to impose 
3(a) 



' ^ ( 5 (2) = 1, a = 3,..., 5), 5 (a) 



p5 + l 1 T 1 1 CK 

1 a 1 



1 5+1 



^ (s(5+l) = l, a = 5 + 2,.., A). (24) 



pal pi 2 \H\"J " ")•••!">/! yv."-V ™ 1 r U|i 

1 1 a 1 2 1 Mo 1 s+1 1 

As in case A, the algebraic relations among the matrices A^ 1 ) (a — 2, N) are fixed by the eigenvalue equation for 
the Hamiltonian ^ when applied to the components of the eigenfunction Q where we have two spins configurations 
| a) (a = 2,...,N) at nearest neighbour sites, but now we have two situations. When we have spins configuration 
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a, (3 < 8 or a, (3 > 8 located at closest positions, and when we have a < 8 ((3 < 8) and (3 > 8 (a > 8). In this last 
case, in order to satisfy the equation Q we need to impose the following constraints: 



(B.i) r^ = r^ = o, r^ = o, t afn 



(b.2) r^-r^-o, rf^r^f - r^r"^, t a01 - tp al - t wi - o, 

where a < <5 (/3 < 5) and (3 > 8 (a > 8), we also obtain the structural constants: 



(25) 



(B.l) S£$(kf\*M) 



S Pa (k i - l3) k {a) ) g(a) v i 2 , v s+ii i(k\ a) +kl f)} ) , (ra i , v i p r n r a /3s 

^aH ' K J > yK yi 21 +i 15+l e 3 ' + l i al + i l / 3~ i ll~ i a/3<' 



(26) 

where we have written for simplicity only the case a < 6 and [3 > 5. On the other hand, when a, /3 < 8 or a, /? > 5, 
we have a similar situation as that considered in the case A and we obtain algebraic relations with Sg t ^ a ,(kj , kj ) 



as in ( 12 ) and coupling constants given by (25 1. The only difference occur when we have a — (3 > 8 where 



(B.2) S^(kf\kf ) ) = l. (27) 

In the sector N — 3 the model B.l generalises the stochastic problem of fully asymmetric diffusion of two kinds of 
particles, whose exact integrability was obtained in |19j through the dynamical matrix product ansatz. The case B.2 
is a new type of integrable model unknown up to now. It is important to mention the B.2 type models can be made 
hermitian and it is related to the A. 4 model by an interchange of labels 1 and a, however since its exactly solutions 
describe different sectors of spin configurations, they are different physical system. 

Finally, as in case A, the momenta fcj (j — 1, n a ) and k^ (j = 1, rib) are fixed by the cyclic property of the 
trace in |5]). We have for both B.l and B.2: 

e -ikfL = A M {k ( a ) Ak (a) }) and e -ikfL = A(»d(jfeC« {jfeW}), (28) 

where A < ' na \kj a \ {k^}) and A^ lb ^(kj b \ {k^}) are the eigenvalue of the transfer matrix defined in (18 1 with 4>(a") — 
9(*'0 L nZiS$%(kf\kM) and m) =9{^) L Y[ZxS%Z{kf\k\ a) ), respectively, and where s£,£„(jfef\*f'>) 



is given by (13 1, (26 1 and (27 1 with rt25 
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DISCUSSION AND CONCLUSIONS 

We solve through a MPA the most general iV-state spin chain with U(1) N symmetry and nearest neighbour inter- 
action in Lunin-Maldacena backgrounds. According to this ansatz, the amplitudes of the eigenfunctions are given in 
terms of a product of matrices where the matrices obey appropriated algebraic relations. This MPA is constructed 
by making a one-to-one correspondence between configurations of spins (or scalar fields in the SYM) and products of 
abstract matrices. Although in Lunin-Maldacena backgrounds all eigenstates have zero momentum, we can generalise 
the MPA including all momentum states [27] • We analyses the Yang-Baxter equation in the N — 3 sector and the 
consistence of the algebraic relations among the matrices defining the MPA and find a new class of exactly integrable 
model unknown up to now (model B.2). This new model is a consequence of the distinguibility of states configurations 
due to the U(l) N symmetry. The study of this new model can be of interest in the context of both AdS/CFT and 
stochastic process. MPA. On the other hand, it will be interesting to sec if it is possible to formulate a MPA for 
the most general one-loop dilatation operator related with an integrable spin model with long-range couplings [39 . 
Another quite interesting problem for the future concerns the formulation of the MPA for the case where we have 
open boundary conditions, as well as for quantum chains with no global conservation law such as the XYZ model, the 
8-vertex model or the case where the quantum chains are defined on open lattices with non-diagonal boundary fields. 
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